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SUMMARY 

Correlation  detectors  operate  by  exploiting  the  similarity  between  two 
waveforms,  one  of  which  is  the  signal  to  be  detected  and  the  other  of  which  is  a 
reference  waveform.  In  practice  the  similarity  occurs  because  both  waveforms 
will  have  been  derived  from  the  same  source.  As  a measure  of  the  similarity 
between  two  signals  we  introduce  here  a 'similarity  function'  which  is  defined 
in  terms  of  the  Fourier  transforms  of  the  signals  rather  than  the  signals  them- 
selves. An  indication  is  given  of  the  limited  validity  of  this  function  as  a 
measure  of  similarity  when  the  signals  are  of  finite  duration  (as  of  course  they 
will  be  in  an  actual  system). 

An  application  of  the  similarity  function  is  discussed  in  which  a linear 
FM  sonar  echo  is  bandpass-filtered  prior  to  being  correlated  with  a replica  of 
the  transmitted  signal.  An  analytical  solution  for  the  similarity  function  in 
terms  of  the  Q-factor  of  the  filter  is  derived  from  which  it  can  be  seen  that  the 
filter  has  little  effect  on  the  detection  capability  of  the  system  provided  that 
the  Q of  the  filter  is  less  than  the  Q-factor  of  the  linear  FM  pulse. 
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INTRODUCTION 


Correlation  detectors  find  application  in  situations  where  the  signal  to 
be  detected  (the  echo)  is  derived  in  some  way  from  some  other  signal,  a replica 
of  which  is  available  as  a reference.  This  causal  connection  between  the  replica 
and  the  echo  suggests  that  the  two  waveforms  should  show  some  similar  features. 

It  does  not  mean  that  they  are  necessarily  identical  in  form.  A measure  of  the 
degree  of  similarity  existing  between  the  two  signals  is  the  entity  which  should 
be  used  as  the  detection  criterion  in  correlation  detection  systems. 


In  section  2 of  this  paper  our  purpose  is  to  examine  the  concept  of 


similarity  between  two  functions  of  time.  Mathematically  the  approach  to  the 

1 . 


problem  is  closely  related  to  that  taken  by  Gabor  in  his  "Theory  of  communica- 
tion". This  is  perhaps  not  surprising  since  in  both  cases  the  aim  is  to  make 
the  best  use  of  the  information  contained  within  a signal. 


1 


In  the  present  case  a 'similarity  function'  is  defined  which,  unlike  the 

cross-correlation  function,  provides  in  general  an  optimum  measure  of  similarity 

between  two  signals.  This  similarity  function  and  the  ambiguity  function  dis- 
2 

cussed  by  Woodward  and  subsequently  used  by  various  other  authors  including 
3 4 

Kramer  and  Harris  and  Kramer  , are  in  fact  the  same  in  the  special  case  where 
the  echo  is  a Doppler-shifted  version  of  the  replica.  The  similarity  function, 
however,  has  wider  application  than  Woodward's  ambiguity  function  and  may  be  used 
quite  generally,  whatever  the  nature  of  the  causal  connection  between  the 
signals . 

An  example  of  how  the  similarity  function  may  be  used  is  discussed  in 
section  3.  The  problem  considered  is  encountered  in  sonar  correlation-detection 
systems  where  the  echo  is  bandpass-filtered  prior  to  being  correlated  with  a 
replica  of  the  transmitted  pulse.  The  effect  of  the  filtering  is  to  diminish 
the  'similarity'  between  the  replica  and  the  echo;  that  is,  the  spectral  compo- 
nents of  the  two  signals  are  no  longer  the  same.  We  use  the  similarity  function 
to  determine  as  a function  of  the  filter  parameters,  by  how  much  the  detection 
capability  of  the  system  is  reduced  as  a result  of  the  filtering  process. 

2 THE  SIMILARITY  FUNCTION 

Let  us  consider  two  real  functions  of  time,  x^(t)  and  x2(t)  , of 
durations  T^  and  Tj  . That  is,  the  amplitude  of  x^(t)  takes  non-zero  values 
only  within  an  interval  of  duration  T^  , beyond  which  x^(t)  is  vanishingly 
small;  and  similarly,  x2(t)  takes  finite  values  only  in  an  interval  of 


U 4^ 


duration  . If  the  leading  edges  of  the  two  functions  occur  at  time  t ■ 0 , 
then 


x^t)  - [u(t)  - u(t  - Tj )] f j ( t ) (1) 

x2(t)  - [u(t)  - u(t  - T2)]f2(t)  (2) 

where  u(t)  is  the  unit  step  function.  And  on  delaying  x2(t)  with  respect  to 
x^(t)  by  a time  t we  obtain  the  signal 

x2(t  - x)  = [u(t  - t)  - u(t  - T2  - x)]f2(t  - t)  . (3) 

If  x > or  t < -T2  there  is  no  overlap  between  x^(t)  and  x2(t  - t). 
Otherwise,  when  the  functions  do  overlap,  the  upper  and  the  lower  bounds,  t2 
and  tj  , of  the  overlap  interval  depend  on  T , T2  and  t in  a way  which  may  be 
determined  by  inspection  of  (1)  and  (3).  The  portions  of  the  two  pulses  within 
the  overlap  region  are  given  by  the  expressions 

x^(t,x)  = [u(t  - tj ) - u(t  - t2>]fj(t)  (4) 

x2(t,x)  = [u(t  - tj^)  - u(t  - t2)]f2<t  - x)  (5) 


and  only  when  x = 0 and  = T2  will  the  overlap  functions  in  (4)  and  (5)  be 
the  same  as  the  signals  in  (1)  and  (2). 


RN 

52 


7 


U 


which  are  just  the  transforms  of  the  two  signals  x^(t)  and  x2<t)  in  the 


absence  of  any  time  delay.  And  the  Fourier  transforms  of  the  two  overlap 
functions  in  (4)  and  (5)  are 


„ 2 

Xj(ju>,i)  - j x'(t,T)exp  - jwtdt  - J f^Oexp  - jwtdt 


(8) 


and 


X2(  ju>,x) 


/ x2(t,x)exp  - jwtdt 

J 


! 


f_(t  - x)exp  - jwtdt  . (9) 


The  transforms  defined  in  equations  (6)  to  (9)  may  be  used  to  construct  the 
function 


J X^( jw,x)X2*( jw,x)dw 


00 

oo 

J |xi(  jo>)  | 2du> 

I |x2(jai)  |2du> 

_0 

0 

(10) 


where  the  asterisk  denotes  complex  conjugate  and  the  integrals,  it  should  be 
noted,  are  taken  only  over  positive  frequencies.  The  implications  of  excluding 
negative  frequencies  from  the  range  of  integration  are  discussed  below. 

Most  practical  correlation-detection  systems  will  be  designed  so  that  over- 
lap losses  are  negligibly  small.  In  other  words,  the  durations  of  the  two 
pulses  will  be  essentially  the  same,  equal  to  T , say;  and  the  value  of  x at 
which  a detection  peak  is  observed  will  be  very  much  less  than  T . Under  these 


conditions,  for  all  frequencies  of  interest,  the  transforms  of  x^(t)  and 


xj(t,x)  are  essentially  the  same;  that  is, 

Xj(ju>,x)  541  X1<ju)  , 

and  the  transform  of  x2(t,x)  is  given  by  the  approximate  relationship 

C ju>,  t ) =“  X2(ju)exp  - ju>x 


(11) 


(12) 


9 


oo 

/ zi 


( ju)Z*( jw)exp  jurrdiu 


*12(t) 


00  T 

1 I ( ju>)  I 2du>  J I Z2(  jco)  1 2du> 


(17) 


where  the  limits  of  integration  now  extend  from  +®  to  . Since  Parseval's 
theorem  states  that 


[ g1(t)g*(t)dt 


OO 

J G^(  jui)G*(  jw)df 
-00 


(18) 


where  g^(t)  and  G^(jw)  are  a Fourier  transform  pair,  we  see  that  an  alterna- 
tive formulation  to  (17)  is 


OD 

/ Z1 


(t)z2(t  - T)dt 


«l2(t) 


tT 

OO  OO  * 

I |zj(t) 1 2dt  I |z  2( t) 1 2dt 


(19) 


The  real  part  of  this  expression  is 


Re  ^^(t) 


00 

/ X1 


(t)x2(t  - T )dt 


r oo  oo 

j x^(t)dt 

j x2(t)dt 

-00 

-00 

TT 


(20) 


which  is  just  the  normalized  cross-correlation  function  between  x^(t)  and 
x2(t)  ; and  the  imaginary  part  is 


10 
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oo 

I |y1(t)x2(t  - t)  - x1(t)y?(t  - t) jdt 


Im  <t^2 ( T ) 


00  oo 

I x*(t)dt  I 


nv 


x2(t)dt 


(21) 


It  is  interesting  to  note  that,  had  we  started  our  analysis  in  the  time  domain, 
we  should  have  chosen  the  correct  function  for  the  real  part  of  ♦ ^(t)  * but 
is  unlikely  that  we  should  ever  have  selected  the  expression  in  (21)  as  the 
imaginary  part  of  ^^(t)  • 

By  the  Schwarz  inequality^  the  following  is  true: 


0 < I *12  ^ T ) I ^ 1 

0 < | Re  *12(t)|  < 1 
0 < | Im  2 ( t ) I ^ 1 


(22) 


and  since  in  general  the  imaginary  part  of  the  similarity  function  is  non-zero, 
it  is  also  true  that 


|*12(t)|  > | Re  4>12(t)| 


(23) 


The  inequality  in  (23)  indicates  why  the  magnitude  of  similarity  function  rather 
than  its  real  part  (the  cross-correlation  function)  should  be  used  as  the 
detection  criterion.  The  difference  between  the  two  functions  is  essentially  a 
difference  in  the  degree  of  sensitivity  each  shows  to  phase  information.  This 
may  be  illustrated  by  considering  as  an  example  the  case  where  the  positive- 
frequency  spectral  components  in  x2(t)  are  all  in  quadrature  with  their 
counterparts  in  x^(t)  and  the  relative  amplitudes  of  the  spectral  components 
are  all  frequency  independent.  The  relationship  between  the  transforms  of  x^(t) 
and  x2(t)  is  then 


X2(jo>)  * 

jXj(  ju>) 

for 

all 

o 

A 

3 

N 

3 

CM 

X 

- jX1(ju) 

for 

all 

O 

V 

3 

(24) 


i 
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and  hence 

*12(0)  - - j ; |*12(o)|  = 1 . (25) 

On  the  other  hand,  we  have 


Re  4^2(0  = 0 for  all  T . (26) 

Clearly,  the  similarity  function  shows  more  tolerance  towards  phase  differences 
between  the  spectral  components  of  x^(t)  and  X2(t)  than  does  the  cross- 
correlation function. 

The  similarity  function  takes  a particularly  simple  form  when  x^(t)  and 
X2(t)  are  the  input  and  output  respectively  of  a linear  system.  If  the  system 
function  is  H(jw)  then  the  relationship  between  the  Fourier  transforms  of  the 
two  signals  is 

X (jw)  = H(jw)X  (jw)  (27) 

and  hence  from  (13) 


CD 

r * 2 

J H ( jw)  |X1(  jw)  | ‘ "exp  jwrdw 


VT) 


TF  • 


J |x^(jw)|2dw  J |H(jw)|2|X1(ju))|2dw 


LO 


(28) 


Assuming  that  the  system  function  is  known,  it  is  necessary  only  to  know  in 
addition  the  power  spectrum  of  the  input  signal  in  order  to  be  able  to  evaluate 
this  expression.  In  the  special  case  where  the  system  function  is  independent 
of  frequency  and  T is  zero,  not  even  the  power  spectrum  is  required,  since 
the  similarity  function  is  simply  equal  to  H*( jw) / | H( jw) | and  | d> ^ ^ ^ 0 ^ I = 1 • 
We  shall  make  use  of  the  formulation  of  the  similarity  function  given  by  (28) 
in  the  following  section. 

3 A PROBLEM  IN  UNDERWATER  ACOUSTIC  SIGNAL  DETECTION 

In  modern  sonar  systems,  linearly-swept  FM  signals  are  used  for  detecting 
objects  whose  radial  velocity  in  relation  to  the  detector  is  too  low  to  shift 
(due  to  the  Doppler  effect)  the  spectral  components  of  the  echo  away  from  the 


frequency  region  in  which  spurious  signals  caused  by  reverberation  predominate. 
Compared  with  its  radar  counterpart,  the  sonar  FM  signal  has  a high  time-bandwidth 
product,  usually  taking  a value  in  the  region  of  several  hundred.  The  implica- 
tions of  this  as  far  as  the  detection  of  a Doppler-shifted  echo  is  concerned  have 

4 

been  discussed  by  various  authors  but  most  notably  by  Kramer  whose  analysis  lays 
the  foundation  of  the  subject.  The  problem  we  discuss  below  is  also  concerned 
with  the  detection  of  a distorted  signal,  but  in  this  case  the  echo  is  a bandpass- 
filtered  version  of  the  linear  FM  replica.  The  discussion  is  of  more  than 
academic  interest  since  in  actual  sonar  systems  it  is  likely  that  the  echo  will 
indeed  be  filtered  in  order  to  remove  the  spectral  components  of  the  noise  lying 
outside  the  frequency  band  of  interest. 

In  order  to  avoid  significant  overlap  losses  the  relaxation  time  of  the 
filter  must  be  very  much  less  than  the  duration  of  the  transmitted  pulse.  For 
most  correlation  sonars  this  condition  would  be  satisfied  if  the  Q-factor  of  the 
filter  were  less  than  about  100.  As  the  Q of  the  filter  vill  be  very  much  less 
than  this  in  practice  (usually  in  the  region  of  20)  in  order  to  maintain  an 
adequately  high  detection  peak,  we  are  wholly  justified  in  assuming  that  the 
durations  of  the  replica  and  of  the  bandpass-filtered  echo  are  effectively  the 
same.  Provided  that  the  time  delay  between  the  replica  and  the  echo  is  very  much 
less  than  the  duration  of  the  two  signals,  the  form  of  the  similarity  function 
given  in  (28)  is  therefore  a valid  measure  for  establishing  the  presence  of  the 
echo. 

We  make  the  further  assumptions  in  the  analysis  that  the  echo  has  been 
reflected  from  a stationary  point  object;  and  that  the  centre  frequencies  of 
the  FM  pulse  and  the  pass-band  of  the  filter  are  the  same.  The  latter  assumption 
is  not  essential  to  the  argument,  but  since  it  leads  to  a simplification  of  the 
algebra  we  use  it  here  in  order  to  illustrate  the  method. 

Let  us  suppose  that  the  transmitted  linear  FM  pulse  has  a centre 
frequency  ojg/2ff  , a frequency-sweep-rate  8 and  duration  T . If  the  function 
x^(t)  of  the  previous  section  is  identified  as  a replica  of  this  pulse,  then 

x^(t)  = [u(t)  - u(t  - T)]f^(t)  (29) 

fx(t)  = sin  j<Kt)f  ; = w0jl  + g(t  - fj| 


where 


(30) 
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From  these  expressions  the  frequencysweep-rate  of  the  pulse  is 


13 


10  “ 0) 

2 1 


where  and  (o^/2tt  are  the  highest  and  lowest  frequencies  in  the  pulse, 

and  Aco/2tt  is  the  pulse  bandwidth. 

In  order  to  evaluate  (28),  the  power  spectrum  of  x^(t)  is  required. 

Fig  1 shows  an  example  of  a computed  spectrum  for  frequencies  close  to  the 
centre  frequency  of  the  replica.  Perhaps  the  most  striking  feature  of  this 
curve  is  its  essentially  rectangular  character.  For  our  present  purpose  we 
ignore  the  structure  on  the  curve  and  approximate  the  spectrum  around  <o  = co^ 
by  a truly  rectangular  distribution.  That  is,  we  put 

S(io)  = I lx1(j“)|2  “ [u(w  " “j)  * u(u>  - «2>]  (32) 

where  S(w)  is  the  power  spectrum  of  the  FM  pulse,  and  the  multiplying  constant, 
tt/Aio  , on  the  right  of  (32)  appears  in  order  to  maintain  consistency  with 
equations  (29)  and  (30). 

The  approximation  in  (32)  should  be  adequate  for  evaluating  the  integrals 
in  the  similarity  function.  For,  if  the  structure  on  the  spectral  density  curve 
were  not  neglected,  its  effect  on  these  integrals  would  more  or  less  average  out 
since  a large  number  of  oscillations  occur  within  the  range  of  integration  (that 
is,  over  the  bandwidth  of  the  pulse). 

By  (32)  the  similarity  function  in  (28)  is  now 


(jw)exp  jtoxdw 


/ l«] 


to)  I ^dio 
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For  a symmetrica  1 response  about 
to  be 

H(  jui) 


a)  ■ u)q  we  choose  the  system  function  in  (33) 


1 + 2jQ 


1 


(a)  * wQ) 


(34) 


The  term  on  the  right  of  (34)  is  the  Lorentz  function.  Strictly  a second 
Lorentz  function  should  be  included  in  (34)  to  ensure  that  H(ju>)  shows  conjugate 
symmetry.  However,  for  frequencies  in  the  region  of  w = , this  second  term 

is  negligibly  small  and  so,  for  our  purpose,  (34)  is  acceptable. 

The  real  and  imaginary  parts  of  the  system  function  are 


Yp) 


1 + 4Q2p2 


Hx(p)  = 


? 2 

1 + 4Q  p 


(35) 


where  we  have  introduced  the  new  variable 


(oj  - 0)Q) 


(36) 


On  recognizing  that  H^tp)  ar*d  Hj(p)  are,  respectively,  even  and  odd  functions 
of  p , we  see  that  | <1>X 2 ^ T ^ I from  (33)  reduces  to  the  form 

1 

2<L 


Qp  I |hr(p)  cos  VT  + Hi^p^  sin  VT}dp 


|*12(t)| 


(37) 


where  = Wq/4w  is  the  'Q-factor'  of  the  pulse.  In  order  to  estimate  the 
magnitude  of  the  principal  maximum  in  this  expression  we  put  t = 0 . In  fact 
this  is  not  an  optimum-  condition  when  the  Q of  the  filter  is  larger  than  Q 

. . P 

The  result  we  obtain  is  therefore  a 'worst  case'  estimate  of  system  performance 

which  will  not  necessarily  be  observed  in  practice. 
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w 


For  zero  time  delay  between  the  replica  and  the  echo,  (37)  becomes 

li 


l*12(0)| 


1 

2Q„ 


f d£_ 

, O + 4( 


2 2 
4Q P ) 


2Q. 


- 

A -1  / Q \ 

u 

7 [%)  ■ 

(38) 

This  expression  is  plotted  in  Fig  2 as  a function  of  Q/Qp  • From  the  shape  of 
this  function  it  is  clear  that  in  practice  it  is  desirable  to  choose  the  Q-factor 
of  the  filter  to  be  no  greater  than  the  Q of  the  pulse.  This  is  not  an 
entirely  unexpected  result  since  it  is  when  Q becomes  larger  than  Q that 

the  spectral  components  of  the  echo  in  the  regions  of  id  = and  id  = 

start  to  suffer  significant  phase  changes  as  a result  of  the  filtering. 

The  analytical  method  described  above  is  valid  for  any  system  function 
provided  that  the  system  to  which  it  relates  is  linear  and,  as  already  stated, 
that  any  relaxation  time  associated  with  the  system  is  short  compared  with  the 
duration  of  the  replica.  For  the  large  class  of  systems  which  can  be  described 
by  a linear  differential  equation  with  constant  coefficients,  H(s)  is  a 
rational  function  of  s ; that  is,  H(s)  is  expressible  as  the  ratio  of  two 
polynomials  in  s . In  most  cases  of  interest  the  order  of  the  denominator  will 
be  higher  than  that  of  the  numerator  (the  system  function  given  in  equation 
(34)  falls  into  this  category);  in  other  words  the  output  of  the  system  will  fall 
to  zero  as  the  frequency  approaches  infinity.  If  H(s)  is  such  a function,  then, 
for  t = 0 at  least,  the  two  integrals  in  (28)  involving  the  system  function 
can  always  be  evaluated  analytically  (assuming  that  (32)  is  used  for  approximat- 
ing the  power  spectrum  of  the  replica),  and  in  general  the  result  for  each 
will  be  the  sum  of  a number  of  logarithmic  functions,  a number  of  inverse  tan- 
gents and  a number  of  rational  functions  of  the  type  A/{(A  - l)(a>  - p)*  } , 
where  A is  an  integer  greater  than  1. 
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